Helix Formation and Folding in an Artificial Peptide 
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We study the relation between a-helix formation and folding for a simple artificial peptide, Alaio- 
Glys-Alaio. Our data rely on multicanonical Monte Carlo simulations where the interactions among 
' all atoms are taken into account. The free-energy landscape of the peptide is evaluated for various 

temperatures. Our data indicate that folding of this peptide is a two-step process: in a first step 
two a-helices are formed which afterwards re-arrange themselves into a U-like structure. 
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I. INTRODUCTION 

The mechanism by which a large class of proteins folds spontaneously into a unique globular shape (!]] has remained 
elusive. Significant new insight was gained over the last few years from the studies of minimal protein models. For 
instance, energy landscape theory and funnel concept proved to be powerful tools for description of the general 
characteristics of folding not only in minimalistic protein models but also for real proteins However, many 

questions on the details of the folding process remain to be solved. For instance, folding of proteins involves one or 
more transitions between different thermodynamic states. The role of these transitions in the folding process is an 
active area of research. An important example for these transitions is the formation of secondary structure elements. 
For the case of a-helices it is long known that there is a sharp transition towards a random coil state when the 
temperature is increased. The characteristics of this so-called helix-coil transition have been studied extensively [gj, 
most recently in Refs. In this paper, we research the relation between helix-coil transition and folding. 

I/*-) ' For this purpose, we have studied an artificial peptide, Alaio-Glys-Alaio, in a detailed representation where the 
I/"") , interactions between all atoms are taken into account. Multicanonical simulations |l3| with large statistics are used 
to evaluate the free energy landscape of our peptide at different temperatures. The encountered transitions are 
further investigated by partition zeros analysis which allows to characterize "phase transitions" in small systems [[l4| . 
Quantities such as the energy, specific heat, helicity and susceptibility were calculated as function of temperature. 
We have neglected in the simulations the interaction of our artifical peptide with the surrounding solvent. While 
this is certainly a crude approximation, it allows us not only to relate our results to our previous studies on helix- 
coil transition in poly-alanine that also relied on gas-phase simulations |c^-[lO|l, but also to study the extend to that 
secondary structure formation and folding are determined by intrinsic properties of the peptide. Our data suggest 
that the peptide in gas-phase folds in a two-step process: in a first step two a-helices are formed in what amounts to 
a first order transition. Afterwards these helices re-arrange themselves into a U-like structure. The second step has 
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q , the characteristics of a second order transition. 
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II. METHODS 



Our investigation of Alaio-Glys-Alaio is based on a detailed, all-atom representation of that peptide. The interaction 
between the atoms is described by a standard force field, ECEPP/2, (as implemented in the program package 
SMMP 0) and is given by: 

Etot = Ec + Elj + Ehb + E tor , (1) 
E c = Y j ^iSl, (2) 
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^ = >;i^-^J. (3) 




E tor = ^Ui{l±cos{ni X i)). (5) 
I 

Here, (in A) is the distance between the atoms i and j, and xi is the /-th torsion angle. The peptide bond angles 
were set to their common value us = 180°. We do not include explicitly the interaction of the peptide with the solvent 
into our simulations and set the dielectric constant e equal to 2. Since the charges at peptide termini are known to 
reduce helix content Jig] , we removed them by taking a neutral NH 2 - group at the N-terminus and a neutral -COOH 
group at the C-terminus. 

Simulation of detailed protein models where the interaction between all atoms are taken into account are extremely 
difficult. This is because the various competing interactions within the molecule lead to an energy landscape char- 
acterized by a multitude of local minima separated by high energy barriers. Hence, in the low-temperature region, 
canonical Monte Carlo or molecular dynamics simulations will tend to get trapped in one of these minima and the 
simulation will not thermalize within the available CPU time. One example of the new and sophisticated algorithms 
p8[ that allow to overcome this difficulty are generalized- ensemble techniques Jl9| , p0[ , and it is one of these techniques, 
multicanonical sampling JT3[, that we used for our investigations. 

In the multicanonical algorithm [ p"3| conformations with energy E are assigned a weight w mu (E) oc l/n(E). Here, 
n(E) is the density of states. A simulation with this weight will lead to a uniform distribution of energy: 

P mu {E) oc n(E) w mu (E) = const . (6) 

This is because the simulation generates a ID random walk in the energy space, allowing itself to escape from any 
local minimum. Since a large range of energies are sampled, one can use the reweighting techniques to calculate 
thermodynamic quantities over a wide range of temperatures T by 

Jdx A(x) w^jEjx)) 
>T ~ Jdx w- l (E{x)) e-^W ' ( ) 

where x stands for configurations. 

Note that unlike in the case of canonical simulations the weights 

w(E) =n- l {E) =e- s(E *> (8) 

are not a priori known. Instead estimators for these weights have to be determined by an iterative procedure |^2],[l2| . 
In our case we needed 500,000 sweeps for the weight factor calculations. All thermodynamic quantities were then 
estimated from one production run of 8, 000, 000 Monte Carlo sweeps which followed 10,000 sweeps for thermalization. 
Our simulations were started from completely random initial conformations (Hot Start) and one Monte Carlo sweep 
updates every torsion angle of the peptide once. At the end of every 10th sweep we stored the ECEPP/2 energies 
Etot, Ec, Elj, Ehb and E tor of the conformation, the corresponding number tih of helical residues and end-to-end 
distance d e - e . Here, we follow previous work [|2j| and consider a residue as helical if its backbone angle (4>,ip) are 
within the range (-70° ± 20°, -37° ± 20°). 

Using the results of our generalized-ensemble simulation, we explored for various temperatures the free energies 

G(n H ,d e - e ) = -IzbT log P(n H , d e -e) ■ (9) 

Here, P(nn, d e — e ) is the probability to find a peptide conformation with values £, d e - e (at temperature T). We chose 
the normalization so that the lowest value of G(uh, d e - e ) is set to zero for each temperature. 

We finally used that the multicanonical algorithm allows us to calculate estimates for the spectral density: 

n(E) = P mu (E)w^ u (E) . (10) 

We can therefore construct the corresponding partition function for our all-atom model of Alaio-Glys-Alaio from 
these estimates by 
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with (3 the inverse temperature, (3 = l/fc^T. The complex solutions of the partition function determine the critical 
behavior of the model and were also studied by us. 

III. RESULTS AND DISCUSSION 

Our peptide, Alaio-Glys-Alaio, is build up out of two chains of each 10 alanine residues connected by 5 glycine 
residues. In previous work p3[Pj9[| we could show that polyalanine has a pronounced transition between a disordered 
coil phase and an ordered state in which the polymer forms an a-helix. For this reason, we expect formation of a-helices 
in our peptide, and the average number of helical residues < n# > is therefore one of the quantities that we have 
measured. < > is displayed in Fig. 1 as a function of temperature, and we observe in this plot two temperature 
regions. At high temperature, few residues are found with backbone dihedral angles ((f>,ip) typical for an a-helix. On 
the other hand, at low temperatures we observe helix-formation, and almost all of the alanine residues are part of an 
a-helix, i.e. have backbone dihedral angles ((f), ip) in the range (—70° ± 20°, —37° ± 20°). The transition between the 
two temperature regions is sharp indicating the existence of a helix-coil transition. The transition temperature The 
can be determined from the corresponding peak in the susceptibility 

X {T)=<n 2 H (T)>-<n H {T)>\ (12) 

which is ploted in the inset of Fig. 1, and we find the transition temperature The = 485 ± 5 K. 

In previous work |2^,|l^l we could show that in polyalanine the formation of a-helices is related to a gain in potential 
energy. For this reason, we display in Fig. 2 the average total ECEPP/2 energy < E to t > and the thermodynamic 
averages of partial energies < Eq >, < Elj >, < Ehb > and < E tor > of Alaio-Gh/5-Alaio as a function of 
temperature. As expected, we observe around our transition temperature The a sharp decrease in the potential energy 
< E to t > that is due to a corresponding decrease in the Coulomb energy < Ec >, Lcnnard- Jones energy < Elj > 
and hydrogen-bond energy < Ehb > ■ The change in < E to t > with temperature can be described by the specific heat 

C(T) = f? - E * ot > - < Etot >2 , (13) 
25 

which we display in Fig. 3. A pronounced change in energy with temperature corresponds to a peak in the specific heat. 
As one can see from Fig. 3, we observe indeed a pronounced peak in the specific heat at a temperature T = 480 ± 10 
K that is consistent with T = 485 ± 5 K, the temperature where the peak in the susceptibility is located. Combining 
both values we obtain as our final estimate for the helix-coil transition temperature The = 483 ± 8. However, we find 
in Fig. 3 also a second, smaller peak in the specific heat at the lower temperature Tf = 265 ± 7 K indicating yet 
another transition. 

The sharp decrease in potential energy corresponding to this second peak is clearly visible in Fig. 2 for the total 
energy < E tot >, however, not for all of the partial energies. Only the Lennard- Jones term < E^j > exhibits also 
a signal for the second transition at Tf. This energy term depends strongly on the overall size of the molecule and 
the change in this quantity indicates a transition between extended and compact structures. Hence, we conjecture 
that the second peak in specific heat at the lower temperature Tf is related to a transition between extended and 
compact structures. A possible measure for such a change is the average end-to-end distance < d e _ e >t- We define 
here < d e _ e > as the distance between N of Alai and O of Ala25, and plot this quantity in Fig. 4. We observe that 
this quantity decreases with decreasing temperature. Below the helix-coil transition The the decrease slows down and 
the curve becomes almost flat at a value of < d e _ e >« 10 A indicating that there is little further change in the 
compactness of the molecule. However, at temperature Tf the end-to-end distance decreases again sharply towards a 
new value < d e - e >= 6.1 A. Hence, Tf marks the folding of the molecule into a defined compact structure with the 
two terminal ends of the peptide close together. 

The transition between extended and a more compact structure can also be seen when we display the free energy 
landscape of our peptide as a function of helicity and end-to-end distance d e _ e . At the temperature T = 480 K 
(which is essentially the helix-coil transition temperature T/ lc = 483 ±8 K) the free energy landscape (displayed in 
Fig. 5a) is flat over a large range of values of uh and d e - e . The 3fcgT contour line surrounds a region where the helicity 
can take values between < nn ^ 20 and the end-to-end distance values between 3 < d e - e < 40, allowing the system 
to move freely between extended and compact configurations, and between helical and coil configurations. On the 



other hand, at the second and lower temperature T = 270 K (which is essentially the folding transition temperature 
Tf = 265 ±7 K) the free energy grows rapidly with decreasing hclicity rin favoring configuration in a small strip 
with 15 < riH ^ 25. Hence, the plot of the free energy landscape in Fig. 5b is limited to values 12 < nu iS 25 of the 
helicity. Here, two regions of minimal free energy can be seen (marked by the 3 ksT contour lines). The first minima 
is found at values of d e - e between 5 and 10 A and 15 < uh < 20 characterising compact structures. A second region 
with slightly lower free energy (see the 1 ksT contour line) is found at much larger values of d e - e between 35 and 40 
A and 20 < tlh < 25 indicating a long stretched a-helix. Both local free energy minima are separated by free energy 
barriers of « 8 ksT that can be overcome by thermal fluctuations. On the other hand, configurations with helicity 
7in < 10 are suppressed by free energy differences of more than 30 kgT. 

Examples for the structures corresponding to the two free energy minima are plotted in Fig. 6. The first one, 
displayed in Fig. 6a is the configuration with lowest energy ever found in our multicanonical simulation of 8,000,000 
sweeps and corresponds to the region in the free-energy landscape at values of <i e _ e between 5 and 10 A and 15 < 
riH < 20. This conformation ('A') consists out of two helixes (made up out of the alanine residues) connected by 
a turn (build out of the flexible glycine residues) towards a U-turn-like structure that is consistent with the small 
value of the end-to-end distance d e - e observed in Fig. 4 for temperatures below Tf. For reference we show in Fig. 6b 
also the configuration ('B') where all 25 residues are part of an a-helix and which corresponds to the second local 
free-energy minimum in Fig. 5b at values of d e - e between 35 and 40 A and 20 < uh < 25 . The dihedral angles 
of both configurations are listed in Table 1. Fig. 7 displays the frequency of both configurations as a function of 
temperature. For T > The neither configuration 'A' nor 'B' are observed. Below that temperature both structures 
appear with similar probability as long as the temperature is higher than Tf. At T = Tf the probability to find 
the maximal helical structure 'B' is with w 30 % highest and decreases after that with decreasing temperature. On 
the other hand, the frequency for the U-turn structure 'A' continues to grow with decreasing temperature. This 
different behavior is due to the energy differences between both structures. The minimal energy conformation 'A' 
has with E to t = —34.7 Kcal/mol a 10.8 Kcal/mol lower potential energy than the extended helix conformation 'B' 
(Etot — —23.9 Kcal/mol). This difference is mainly due to the Lennard-Jones terms: Elj = —132.5 Kcal/mol for 
'A' vs. Elj = —118.9 Kcal/mol for 'B'. The gain in Elj is in part compensated by the hydrogen-bonding terms: 
Ehb = —30.2 Kcal/mol for 'A' vs. Ehb = —34.7 Kcal/mol for 'B'. Coulomb and torsion energies differed little between 
the two configurations: Ec — 126.7 Kcal/mol and E tor = 1.3 Kcal/mol for conformation 'A' vs. Ec = 126.4 Kcal/mol 
and E tor = 3.3 Kcal/mol for conformation 'B'. 

It is an interesting question whether our two observed transitions (occuring in a finite and small system) can be 
related to phase transitions which in a strict sense are defined only for macroscopic (that is very large) systems. In 
order to study this question we have calculated the complex zeros j3 — > Re(/3) + it of the partition function Z(j3) of 
our molecule. In the case of a temperature driven phase transition, we expect that the complex zeros j3j, (j = 1, 2, ...) 
(or at least the ones close to the real axis) condense for large enough system size on a single line. As the system size 
increases, those zeros will move towards the positive real /3-axis and the corresponding value is for large system size 
the inverse of the physical critical temperature T c . Crucial information on phase transitions can be obtained from 
the way in which the first zero approaches the real /3-axis. However, such an analysis depends on the extrapolation 
towards the infinite large system and does not allow characterization of the situation in small systems such as Alaio- 
Gly5-Alaio . One possible extension of the above ideas to "phase transitions" in biological molecules and other small 



systems is the classification scheme by Borrmann et al. |14|. In this approach one computes the discrete line density 



of zeros as an average of the inverse distances between neighboring zeros, 

^ Tfe ) = \ ( Tfl — \ — i + Vr — > ( 14 ) 

2 \\Pk- Pfe-il IPfc+i - Pk\J 
and approximates </>(r) by a simple power law 4>{t) ~ r Q . Taking the first four complex zeros, one obtains 

ln0(r 3 ) -ln0(r 2 ) 

a= : : . (15) 

m r 3 - in t 2 

With a second parameter 7, related to the crossing angle of this line with the real axis, 

7 =[Re(/3 2 )-Re(/3 1 )]/(r 2 -ri), (16) 

and following the classification scheme by Grossmann and Rosenhauer J24|,|2j| , phase transitions can now be classified 
according to the values of these two parameters: for a < and 7 = one has a phase transition of first order, it is of 
second order if < a < 1 and arbitrary 7, and for a > 1 and arbitrary 7 one has a higher order transition. We have 
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evaluated the usefulness of this approach both for spin systems and polyalanine chains |2qj27|] . Preliminary results 
for Alaio-Gly5-Alaio are also listed in Ref . |27j] ■ 

For Alaio-Gly5-Alaio, we find two lines of complex zeros. The corresponding first four zeros for each characteristic 
line are listed in Table 2. Our error estimate is based on the jackknife method |2£| with 16 bins. These lines lead 
to two critical temperatures The — 480 K and Tf = 271 K (estimated from the real part of (3) that agree with the 
corresponding values The = 483 ±8 K and Tf = 265 ± 7 K found by us above with different methods. 

Using these zeros we have calculated the parameters a and 7 that characterize in the Borrmann et al. approach 
phase transitions in small systems. For the first transition, at T = 480 K, we find a = 1.1(1.5) and 7 = —0.4(2). 
The errors reflect large fluctuations in the values of the two parameters a and 7 that do not allow us to determine 
whether the helix-coil transition is a weak first order or a strong second order phase transition. This problem was 
also observed in our earlier work on polyalanine |9|^^| where we were also not able to establish clearly the order of the 
helix-coil transition. However, our results illustrate the strength of this transition that also leads to the pronounced 
peak in the specific heat observed in Fig. 3, and suggest a nucleation mechanism for a-helix formation. Our data are 
more decisive in the case of the second transition, at T — 265 K, which marks the collapse and folding of the peptide. 
Here we find a = 0.32(8) and 7 = 0.36(2). These values indicate a second-order transition which is consistent with 
what one would expect for a transition between extended and compact structures and imply that collapse and folding 
of the Alaio-Gly5-Alaio is connected with long range correlations between the residues. 

Our above analysis of the thermodynamics of our peptide suggests that Alaio-Glys-Alaio folds in a 2 step process. 
The first step is the formation of a-helices and can be characterized by a helix-coil transition temperature The = 483 ±8 
K. The formation of a-helices then restricts the possible configuration space. Energetically most favorable is the folding 
of two a-helices (made out of the alanine residues) into a hairpin. This second step can be characterized by a lower 
folding temperature Tf — 265 ± 7 K. Note that this folding temperature is in the biological relevant temperature 
regime while helix-formation can also happen at much higher temperatures. The above described two step folding of 
our artificial peptide is reminiscent of the well known framework [^9|J3C| 1 and collision-diffusion model 1 31 of folding 
which also propose that local elements of native local secondary structure form independently of tertiary structure. 
These elements diffuse until they collide and coalesce to give a tertiary structure. In our case, the temperature region 
of 265 — 480 K is the one where the thermal energy of the molecule does not allow coalescing of the helix-fragments 
that therefore form and decay. Some stabilization happens when these fragments try to form one extended helix, 
however, the inherent flexibility of the glycine residues, connecting the two alanine chains, and the gain in Lennard- 
Jones energy lead instead at temperatures below Tf to a U-turn-like bundle of two (antiparallel) a-helices connected 
by a turn of glycine residues as the most stable structure. Note that this picture is consistent with energy landscape 
theory and funnel concept (|||. Fig. 5b depicts the appearance of a folding funnel at T = 270 K towards our "native 
structure" 'A'. The competing structure 'B', that at this temperature has a slightly higher (« lfcgT) free energy ( 
see the 1 fcgT contour line for structure 'A' that is missing for conformer 'B'), acts as a local trap. However, the free 
energy barriers of w 6 fc^T can be overcome at this temperature by thermal fluctuations. Below that temperature 
the relative weight of structure 'B' decreases (see Fig. 7) and its free energy difference to 'A' increases: the energy 
landscape becomes even more funnel like (data not shown) . The energy landscape of Fig. 5b allows for a multitude 
of folding pathways that all, however, follow the above described two-step process. 

An interesting question is how general the above obtained results are. A direct comparison with experimental 
data is difficult since solvent effects were neglected in the simulation of Alaio-Glys-Alaio and most experiments study 
solvated peptides. An exception are the techniques developed by Jarrold and collaborators for examination of gas- 
phase conformations of proteins and peptides |32[ . An experimental study of Alaio-Glys-Alaio using these techniques 
is now planed. The authors are not aware of experimental results for the solvated peptide. In order to compare our 
work with experiments of other short helical peptides |3^JU], we therefore started now simulations of Alaio-Gly5- 
Alaio where the solvation effects are approximated by a solvent accessible surface term Q. This will allows us also 
to test the dependency of our results on the solvation model. Simulating another, slightly more complicated, artifical 
peptide, Alaio-Glys-Alaio-Glys-Alaio that presumably will fold in a three-helix bundle, will allow in addition a direct 
comparison with recent experimental work by Myers and Oas |34| where the relation between helix-formation and 
folding was studied for the 58-residue B domain of protein A. 



IV. CONCLUSION 



In summary, we have performed multicanonical simulations with high statistics of a simple artificial peptide, the 
25 residue Alaio-Glys-Alaio- We found that this peptide folds into a specific structure that is determined solely by 
the intrinsic properties of the molecule (since solvent interactions are absent in our simulations. In gas-phase, the 



peptide exhibits two characteristic transitions. At The = 483 ± 8 K we observed a helix-coil transition that is either a 
weak first order transition or a strong second order transition. Our results indicate that there is a second transition at 
Tf = 265 ± 7 K, the folding transition, that is second order-like. These results suggest that folding of this peptide in 
gas-phase is a two-step process. In a first step, the alanine residues form independently helical segments which then 
afterwards in a second step assemble to a U-turn like structure of two antiparallcl a-helices connected by a turn. By 
using an implicit solvent model in our simulations we started now to investigate whether the final structure or this 
two-step process changes in the present of water. 
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Table I. Dihedral angles for the two configurations shown in Fig. 5 



State 




Conformer A 






Conformer B 




Etot (Kcal/mol) 






-34.7 






-23.9 




Residue 




<t> 




X 






X 


Ala-f 




-101.7 


154.4 


-178.1 


-55.2 


-54.9 


61.5 


Ala-2 




-65.4 


-35.5 


-173.8 


-70.6 


-37.8 


49.2 


Ala-3 




-68.2 


-36.1 


-61.7 


-60.4 


-30.8 


-178.3 


Ala-4 




-70.6 


-39.0 


-69.7 


-78.6 


-37.6 


-43.4 


Ala-5 




-64.3 


-40.6 


-54.1 


-66.5 


-28.8 


176.5 


Ala-6 




-66.3 


-38.3 


175.1 


-70.6 


-50.4 


-55.8 


Ala-7 




-69.5 


-31.9 


-172.6 


-64.2 


-32.6 


-178.6 


Ala-8 




-75.5 


-31.6 


57.2 


-68.6 


-43.4 


-58.3 


Ala-9 




-62.6 


-45.2 


-172.9 


-67.8 


-4L0 


-164.8 


Ala-fO 




-73.4 


-54.4 


-59.0 


-60.8 


-43.3 


53.4 


Gly-ff 




-91.9 


65.0 




-64.3 


-40.5 




Gly-12 




156.2 


-80.6 




-64.9 


-44.9 




Gly-f3 




151.5 


-176.6 
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Table II. 


Partition function 


zeros for the two transitions observed for Alaio 


-Glys-Alaio- 




Re(/3i) 


Tl 


Re(/3 2 ) 


T2 


Re(/?s) r 3 


Re(/3 4 ) 




1.0463(46) 


0.1307(53) 


f.0f44(99) 


0.2ff2(55) 


1.051(44) 0.310(36) 


1.055(46) 


0.347(41) 


1.855(21) 


0.263(12) 


1.991(25) 


0.637(4f) 


2.057(53) 0.923(43) 


2.070(43) 


1.229(55) 
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Figure Captions: 



Fig. 1 The average helicity < nu > (T) of Alaio-Glys-Alaio as a function of temperature T. The corresponding 
values of the susceptibility x(T) are ploted in the inset. Our data rely on a single multicanonical simulation of 
8,000,000 Monte Carlo sweeps. 

Fig. 2 The average ECEPP/2 energy < E to t > and the corresponding partial energy terms, Coulomb energy < Eq >, 
Lennard- Jones term < Elj >, hydrogen-bond energy < Ehb > and torsion energy < Et or >, as a function of 
temperature T for Alaio-Glys-Alaio- 

Fig. 3 The specific heat C{T) of Alai -Gly 5 -Alaio as a function of temperature T. 

Fig. 4 The average end-to-end distance < d e _ e > (T) as a function of temperature T for Alaio-Gly 5 -Ala 10 . 

Fig. 5 Free energy landscape of Alaio-Glys-Alaio as a function of helicity nn and end-to-end distance d e - e at 
temperature (a) The = 480 K and (b) Tf = 270 K. The contour lines are drawn in multiples of fcgT indicated 
in the key. 

Fig. 6 Lowest-energy conformation 'A' of Alaio-Gly5-Alaio (a) and the conformation 'B' (which has the largest 
helicity) (b) as found in our multicanonical simulation of 8,000,000 Monte Carlo sweeps. 

Fig. 7 Relative weight of the lowest-energy conformation 'A' and conformation 'B' (which has maximal helicity) as 
a function of temperature. 
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(Fig.l) 1. The average helicity < nn > (T) of Alaio-Glys-Alaio as a function of temperature T. The corresponding 
values of the susceptibility x(T) are ploted in the inset. Our data rely on a single multicanonical simulation of 8,000,000 
Monte Carlo sweeps. 
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(Fig. 3) 3. The specific heat C(T) of Alaio-Glys-Alaio as a function of temperature T. 
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(Fig. 4) 4. The average end-to-end distance < d e _ e > (T) as a function of temperature T for Alaio-Glys-Alaio. 
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(Fig.5a) 5. Free energy landscape of Alaio-Glys-Alaio as a function of helicity nn and end-to-end distance d e _ e at 
temperature (a) T hc = 480 K and (b) T f = 270 K. 
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(Fig.6a) 7. Lowest-energy conformation 'A' of Alaio-Glys-Alaio (a) and the conformation 'B' (which has the largest 
helicity) (b) as found in our multicanonical simulation of 8,000,000 Monte Carlo sweeps. 
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(Fig.7) 9. Relative weight of the lowest-energy conformation 'A' and conformation 'B' (which has maximal helicity) 
as a function of temperature. 
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